o 

(N 



Guided nucleation of superconductivity on a graded magnetic substrate 
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We demonstrate the controlled spatial nucleation of superconductivity in a thin film deposited 
on periodic arrays of ferromagnetic dots with gradually increasing diameter. The perpendicular 
magnetization of the dots induces vortex-antivortex molecules in the sample, with the number of 
(anti)vortices increasing with magnet size. The resulting gradient of antivortex density between the 
dots predetermines local nucleation of superconductivity in the sample as a function of the applied 
external field and temperature. In addition, the compensation between the applied magnetic field 
and the antivortices results in an unprecedented enhancement of the critical temperature. 

PACS numbers: 73.23.-b, 74.78. Na 
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The enhancement of the superconducting critical pa- 
rameters is arguably the most prominent objective of 
mesoscopic superconductivity. For example, it is already 
well established that critical magnetic fields rise in sub- 
micron samples to values substantially larger than those 
found in bulk materials [1|, and both critical field and 
current can be enhanced by nanoengineered periodic [2| 
and quasi-periodic pinning [3]. 

The best current method of preventing the deteriora- 
tion of the superconducting state under external mag- 
netic fields consists of locally counteracting the applied 
field by an array of ferromagnets with perpendicular 
anisotropy. These field-compensation effects can lead to 
a substantial increase of the Tc{H) boundary in a field 
range determined by the magnetization of the dots [J]. 
More precisely, the maximum of the Tc{H) boundary 
shifts away from zero field, to a value determined by 
the number of antivortices in the sample [5|- Namely, 
if magnets are capable of generating vortex-antivortex 
molecules in the superconductor ;6!] , an applied field that 
provides exactly the same number of vortices will lead 
to a complete vortex-antivortex annihilation between the 
dots, which will maximize Tc{H). 

In an attempt to achieve the ultimate field compensa- 
tion, here we study a thin superconducting film evapo- 
rated on top of a square array of out-of-plane magnetized 




FIG. 1: (color online) (a) A square lattice of magnetic dots 
with linearly increasing size on top of a superconducting film 
and an oxide layer, (b) The stray field of each dot {km) may 
induce vortex-antivortex molecules in the superconductor 



dots o/ variable size, i.e. parallel periodic rows of dots in- 
creasing in diameter [see Fig. [Ua)]. We expect that such 
gradual magnetic texturing will enhance superconductiv- 
ity in different parts of the sample at different applied 
fields, and therefore should lift the Tc{H) characteristics 
as a whole. Since vortex-antivortex (V-Av) physics is 
fundamentally important in such a system where differ- 
ent dots induce different numbers of V-Av pairs [see Fig. 
[Ijb)], we also investigate the resulting structure of V-Av 
complexes. 

Our sample is a 50 nm thick Al film with critical 
temperature Tco = 1.350 K. Experimentally estimated 
values of coherence length and penetration depth are 
^(r=0)=117 nm and A(0)=98 nm, respectively. A 5 nm 
thick Si insulating layer separates the Al film from the 
underlying square array of circular magnetic dots with 
W = 2 fim lattice spacing [see Fig. HJa)], so that any 
proximity effect can be neglected. The ferromagnetic 
dots have diameters of 0.2, 0.4, 0.6, 0.8 and 1.0/im and 
consist of a 2.5 nm Ft buffer layer covered with a [0.4 
nm Co/ 1.0 nm Pt]io multilayer with magnetization per- 
pendicular to the sample surface. Our numerical simu- 
lations are performed both within the Ginzburg-Landau 
(GL) formalism and by molecular dynamics (MD) using 
vortex- vortex and magnetic dot- vortex interaction poten- 
tials derived from London theory JJ . For thin supercon- 
ducting films (always effectively type-II) the results of 
these two theories should converge. 

Vortex-antivortex complexes To analyze theoretically 
the vortex-antivortex configuration and their influence 
on the S/N phase boundary, we employ the GL theory. 
The two GL equations for the order parameter ^ and 
the vector potential A are solved self-consistently, as de- 
tailed in Ref. [6[. The ground-state vortex configura- 
tions are determined by comparing the Gibbs free energy 
j^/Tq = y-i /(2(A - Ao)j - I^Hdr of all found stable 
vortex states, where Jx) = H^/8n, and all quantities are 
dimensionless. Here, j denotes the supercurrent density. 
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FIG. 2; (color online) The GL free energy of the supercon- 
ducting state as a function of the magnetization of the mag- 
netic dots (field-cooled regime), (a-e) Local magnetic field 
distribution in the sample for representative vortex-antivortex 
configurations (N = 1,2,4,6, 13 antivortices per row respec- 
tively), and comparison with MD simulations (f,g). 



Aq the applied vector potential, and V the sample vol- 
ume. In what follows, most of the calculated quantities 
will be expressed (back) in real units for direct compari- 
son with experiment. 

Fig. [2] shows the calculated evolution of the ground- 
state energy of the superconducting film as a function 
of the magnetization of the magnetic dots, for other pa- 
rameters taken as in the experiment (d, I, and dm denote 
the thickness of the film, oxide and the magnetic dots 
respectively). As schematically shown in Fig. [IJb), each 
dot generates an inhomogeneous stray field which favors 
the nucleation of vortex-antivortex (V-Av) pairs. For a 
given magnetization, V-Av pairs first nucleate under the 
largest dots, due to their larger magnetic moment [see 
Fig. [Ua)]. In interaction with other dots, antivortices 
occupy interstitial positions, while vortices remain under 
the largest dots at a somewhat off-center position. With 
increasing magnetization, more V-Av pairs are induced 
in the sample, and they gradually appear around smaller 
dots as well [see Figs. [2Ib-d)]. Increased density of an- 
tivortices between the magnetic dots leads to their 'crys- 
tallization' into an interstitial lattice [Fig. [2le)]. Results 
of the GL theory are nicely corroborated by MD simula- 
tions [see Figs. [2jf,g)]. 

Two important conclusions follow from the above re- 
sults. First, our gradiated magnetic texture leads to a 



gradient of antivortex density along the sample (for vor- 
tices as well, but they are localized under the magnetic 
dots). Second, depending on the spacing between the 
magnetic dots and the number of nucleated antivortices, 
some of them can penetrate even the regions free of nu- 
cleated vortex-antivortex pairs [see area around smallest 
dots in Fig. [UJe)]. This shows that antivortices can de- 
tach from the original magnetic dot and move towards 
other dots, guided by an overall magnetic potential. 

The latter conclusions are even more evident if realistic 
structural defects are introduced in the calculation. In 
Fig. ISlJa) we show the vortex-antivortex state obtained 
for M — 300 kA/m, where we allowed for minor devia- 
tions of the size of the magnets (< 1 nm) compared to the 
experimentally suggested values. As a result, some dots 
with the same nominal size generate different number of 
V-Av pairs in the sample and consequently the overall 
magnetic profile becomes less uniform and antivortices 
penetrate easier the areas around the dots of insufficient 
moment for the nucleation of V-Av pairs. Interestingly, 
the magnetic profile in the sample changes considerably 
with increasing temperature [see Figs. EJb-d)]. Fig. [3] 
shows that higher temperature leads to: (i) appearance 
of new V-Av pairs around some dots; (ii) merging of con- 
fined vortices under the dots into a giant-vortex [8|; (iii) 
increased rigidity of the antivortex lattice between the 
magnetic dots; and (iv) gradual destruction of supercon- 
ductivity, starting from areas around the largest dots. 
Molecular dynamics simulations disregard the suppres- 
sion of superconductivity under the magnetic dots and 
do not take into account the large vortex cores - as a 
consequence they fail to reproduce the GL results at tem- 
peratures above 0.8 Tc. 
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FIG. 3: (color online) The evolution of the Cooper-pair con- 
densate (and its energy) under the gradiated array of mag- 
netic dots as a function of temperature (for Af — 300 kA/m). 



Controlled nucleation of superconductivity The above 
findings indicate that with increasing temperature super- 
conductivity is longest preserved in the region around 
the smallest magnetic dots. However, it is already a 
well-known fact that nanostructuring of superconduct- 
ing films by perpendicularly magnetized dots leads to a 
shift of the H — T phase boundary, a phenomenon known 
as field-induced-superconductivity [J] . This phenomenon 
is asymmetric with applied external field, i.e. fields par- 
allel to the moment of magnetic dots enhance supercon- 
ductivity (due to stray field compensation and vortex- 
antivortex annihilation between the dots) while applied 
antiparallel field only suppresses superconductivity fur- 
ther. As a result, the H — T boundary shifts with re- 
spect to H = Q, and exhibits maximal Tc{H) at a field 
determined by the magnetic moment of the dots. In our 
present sample however different rows of magnetic dots 
have different magnetic moments, and it is by no means 
obvious how the H — T boundary will look like. 

The results of our transport measurements and GL 
simulations for the superconducting/normal phase H — T 
boundary are shown in Fig. 3] In the theoretical simula- 
tions, the nucleation of superconductivity was monitored 
locally, i.e. maximal |^p > 10^^'^ was declared as a 
superconducting state. In the case of negative applied 
field, we obtained an identical phase boundary for both 
magnetized and demagnetized dots (dashed curve in Fig. 
S]). In this case, since superconductivity always nucle- 
ates along the line between the rows of the two small- 
est dots, Tc{H) shows an approximately linear behav- 
ior. Deeper inside the superconducting phase (lower T) 
the applied magnetic field induces additional antivortices 
in the sample, increasing the configurational complexity. 
As an illustration, the inset of Fig. g] at T = 0.86 Tc 
and H = -ff_4, i.e. matching field adding 4 antivortices 
per magnetic dot, shows the resulting structure of both 
antivortex molecules around the magnetic dots and an 
antivortex lattice between them. 

On the other hand, for applied positive magnetic field 
we observe a clear enhancement of superconductivity in 
our sample. As detailed in Ref. [5[, the highest Tc in 
magnetically structured superconducting films is always 
obtained at the applied field matching the number of an- 
tivortices per periodic unit cell of the sample. Since our 
sample contains magnetic dots of five different sizes and 
correspondingly different number of induced antivortices 
around them (iV=l-|-2-f 2-1-5-1-8), we achieve the latter 
condition for maximal Tc at several values of applied 
magnetic field! As a consequence the overall maximum 
of Tc remains at iJ = 0, but the Tc{H) boundary is 
lifted tremendously as compared to the reference sam- 
ple. The superposition of compensation phenomena in 
different parts of the sample leads to several cusps in 
Tc{H) [9|, but the general trend is clear - as the applied 
magnetic field is increased, superconductivity nucleates 
first in the regions of the sample with matching antivor- 



Tb«;urelitAlCiik;u3Htion 

Demagnelized dots 



0.8S 




FIG. 4: (color online) The superconducting/normal (S/N) 
H — T phase boundary obtained by transport measurements, 
compared to the case of demagnetized dots and the theoretical 
data (for M = 450 kA/m, A(0) = 100 nm, ,^(0) = 110 nm). 
The insets show the Cooper-pair density plots at (H,T) values 
denoted by open dots. 



tex density. As illustrated in Fig. |4l the nucleation line 
shifts from smallest (at H = Hi) towards largest dots (for 
H = Hg,)- Such guided nucleation of superconductivity 
is not only a unique fundamental phenomenon, but it is 
of value for potential multi-channel applications, and, as 
shown in Fig. 21 leads to an unprecedented enhancement 
of superconductivity in an applied magnetic field. 

To summarize, thanks to the specific design of the ar- 
rays of magnetic dots with varied size, and the result- 
ing gradient in antivortex density in the superconducting 
film, we achieved the highest overall critical temperature 
of superconducting films in the presence of applied mag- 
netic field. We also show that nanoengineered arrays of 
magnetic dots enable the control of local nucleation of 
superconductivity by simply changing the magnitude of 
applied magnetic field. 
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